Canonical Polygon Queries on the Plane: A New 

Approach 

S. Sioutas, D. Sofotassios, K. Tsichlas, D. Sotiropoulos, P. Vlamos 
Ionian University-Department of Informatics, Computer Technology Institute, Aristotle University of 

Thessaloniki-Informatics Department, Greece 
Email: sioutas@ionio.gr, sofos@cti.gr, tsichlas@csd.auth.gr, {dgs,vlamos}@ionio.gr 



Abstract — The polygon retrieval problem on points is the 
problem of preprocessing a set of n points on the plane, 
so that given a polygon query, the subset of points lying 
inside it can be reported efficiently. It is of great interest in 
areas such as Computer Graphics, CAD applications, Spatial 
Databases and GIS developing tasks. In this paper we study 
the problem of canonical k -vertex polygon queries on the 
plane. A canonical fc-vertex polygon query always meets 
the following specific property: a point retrieval query can 
be transformed into a linear number (with respect to the 
number of vertices) of point retrievals for orthogonal objects 
such as rectangles and triangles (throughout this work we 
call a triangle orthogonal iff two of its edges are axis- 
parallel). We present two new algorithms for this problem. 
The first one requires 0(n log 2 n) space and 0(k t l ° g l0 g n +A) 
query time. A simple modification scheme on first algorithm 
lead us to a second solution, which consumes 0(n 2 ) space 
and O(fc ; '° ; og„ + A) query time, where A denotes the size 
of the answer and k is the number of vertices. The best 
previous solution for the general polygon retrieval problem 
uses 0(n 2 ) space and answers a query in 0(fclogn + A) 
time, where k is the number of vertices. It is also very 
complicated and difficult to be implemented in a standard 
imperative programming language such as C or C++. 

Index Terms — Algorithms and Complexity, Data Struc- 
tures, Computational Geometry, Spatial Databases. 

I. Introduction 

Given a set S of n points on the plane, the problem 
of retrieving a subset S' 6 S that lie in the interior of 
a planar geometric object is of great interest in the areas 
of Computational Geometry, Spatial Databases, Computer 
Graphics, CAD and GIS applications. The efficiency of 
the solutions presented so far depends on the existence 
or not of orthogonality on the query object, which means 
that not all the line segments forming the query figure are 
vertical or horizontal. 

A range tree [20] for a example, is powerful enough 
to support windowing of points (i.e. the query object is 
an arbitrary axis-parallel rectangle) in 0(\ogn + A) time 
using O(nlogn) space. The problem becomes harder as 
the complexity of the query object increases (i.e. triangle, 
quadrilateral, arbitrary polygon), and there is no full 
orthogonality. 

It is important to notice that the point retrieval problem 
for simple polygons is an interesting problem for many 
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application areas. In medicine for example, the term ROI, 
which stands for Region Of Interest, is widely used by 
physisians in order to indicate a polygonal region of 
arbitrary complexity on their scene (i.e. Radiology Im- 
age). Many information systems developed so far support 
retrieval queries at such a region by first computing its 
bounding rectangle. Then, they report all the points inside 
the rectangle and finally these points are filtered so that 
only the points inside the region remain. 

Willard [15], was the first to present a solution for the 
point retrieval problem for simple fc-vertex polygons. It 
uses 0(n) space and the query time is O(n 77 + A). 
Edelsbruner and Welzl [10], reduced the query time 
to O(n - 69 + A). A faster algorithm was presented 
by Edelsbruner, Kirkpatrick and Maurer [9] that uses 
0(fclogn + A) query time and 0(n 7 ) space. Cole and 
Yap [6], presented a method using 0(n 2 / \ogn) space 
and 0(fclog7jloglogn + A) time. Finally, Paterson and 
Yao [13] have presented the best known solution (for a 
simple arbitrary polygon). This solution uses 0(n 2 ) space 
and answers a query in 0(fclogn + A) time. 

In this work we consider the polygon retrieval problem 
on points in special case of canonical polygons which 
always satisfy the following strict property: a retrieval 
query on a set of points can be transformed to a linear 
number (on the number of their vertices) of queries on 
orthogonal objects such as rectangles and triangles. We 
call a triangle orthogonal iff two of its edges are axis- 
parallel. We present two solutions: The first one requires 
0(nlog 2 n) space and 0{k l l °^ lo ^ n + A) query time. The 
second one consumes 0(n log n) space and 0(k ^"f^ + 
A) query time, where A denotes the size of the answer 
and k is the number of vertices. 

The best previous solution for the general polygon 
retrieval problem uses 0(n 2 ) space and answers a query 
in 0(k \ogn + A) time, where k is the number of vertices. 

This paper is organized as follows. In Section [TT] we 
briefly introduce some preliminary data structures. In 
Section [Til] we give the details of our algorithms. In 
Section [TV] we present some special extensions for the 
general polygon retrieval problem based on algorithm 
of Paterson and Yao for which we shortly introduce 
the fundamental notions. Finally, some conclusions and 
further extensions concerning this problem are considered 
in Section [Vl 



II. Preliminary Data Structures 

A. Fusion Trees 

At ACM STOC 1990, Fredman and Willard [11] sur- 
passed the comparison-based lower bounds for sorting and 
searching using the features in a standard imperative pro- 
gramming languages such as C. Their key result was an 
0(logn / loglogn) time bound for deterministic searching 
in linear space. The time bounds for dynamic searching 
include both searching and updates. Since then much 
effort has been spent on finding the inherent complexity 
of fundamental searching problems. 

B. Planar Point Location in Sublogarithmic Time 

Timothy M. Chan [1] and Mihai Patrascou [18] ex- 
tended the 1 -dimensional fusion tree [1 1] to 2-dimensions, 
in order to handle the point location problem in 
0(logn/ loglogn) time and linear 0(n) space. 

C. Half Plane Range Query 

The half plane range query problem is the problem of 
reporting all the points in a set S of n points on the plane 
that lie on a given side of a query line L. This section 
combines the method presented by Chazelle, Guibas, and 
Lee [5] that achieves 0(logn + A) query time and linear 
space (using the notion of duality) with the best current 
method for planar point location problem [1], [18]. The 
main steps of algorithm [5] are the following: 

« Preprocessing 

1) Partition S into a set of convex layers: 

(a) Define Si as the convex hull of all the 
points currently in S 

(b) Remove the vertices of Si from S 

(c) Increment i, repeat the process 

The time cost is 0(n\ogn) while the space 
complexity is 0(n), using a technique that 
computes convex hulls in a dynamic environ- 
ment [2]. 

2) Augment the set of layers building vertical con- 
nections as follows: for each vertex w of layer 
Si, keep a pointer to the two edges immediately 
above and below w. This clearly uses 0(n) 
extra space. 

3) Using duality, the transformation of each vertex 
w into its corresponding line maps each layer 
into another convex polygon. The produced 
mapping is organized into a point location 
structure, occupying 0(n) space. 

> Query Processing 

1) Given a query line L transform it into its 
corresponding dual point P^. 

2) Apply a planar point location algorithm for the 
point Pl in the properly organized structure. 
This determines the innermost layer among the 



layers containing the point Pl- Thus, in the 
dual mapping it determines the innermost layer 
among the layers that L intersects. Call this 
layer neighboring. Using the best current point 
location algorithm, this costs 0(logn/ loglogn) 
time. 

3) Using the pointers mentioned at step [3] of the 
preprocessing procedure, it is easy to report one 
vertex lying at the query half plane for each 
layer which encloses the neighboring one. 

4) Traverse each layer from each of the vertices 
reported across the part of the layer inside the 
half plane. Report the vertices traversed. 

Clearly, these steps lead to an algorithm for answer- 
ing half plane range queries using 0(n) space and 
Oilogn I loglogn + A) query time. We use this method 
in order to answer orthogonal triangle range queries on 
points. 

D. Priority Search Tree 

In this subsection, we briefly review the priority search 
tree of McCreight [16]. Let S be a set of n points on the 
plane. We want to store them in a data structure, so that 
the points that lie inside a semi-infinite strip of the form 
(— oo, 6] x (— oo, c], can be found efficiently. 

The priority search tree is a binary search tree over the 
^-coordinates of the points. The root of the tree contains 
the point p with the minimum y-coordinate. The left (resp. 
right) subtree is recursively defined for the set of points 
in S — {p}. The set S — {p} is partitioned equally into 
the two subtrees of the root. As a result, it is easy to see, 
that a point is stored in a node on the search path from 
the root to the leaf containing its .x-coordinate. 

Queries with ranges that are half-infinite in both x and 
y directions are also known as quadrant range search. 
To answer a quadrant range query, we find the O(logn) 
nodes in the search path Pi, for point b. Let Lf, be the 
left children of these nodes that do not lie on the path 
(see Figure [TJ. In O(logn) time, the points of the nodes 
of Ph [J Lb that lie in the query-range can be determined. 
Then, for each node of L& storing a point inside the range 
query, its two children are visited and checked whether 
their points lie in the range. This procedure continues 
recursively, as long as points in the query-range are found. 

The correctness of the query algorithm is proved as 
follows. First, observe that nodes to the right of the search 
path, have points with ^-coordinate larger than b and 
therefore lie outside the query-range. The points of 
may have x-coordinate larger than b or they may have y- 
coordinate larger than c. In any case, they are not reported. 
The nodes of Lf, and their descendants have points with IE- 
coordinate smaller than 6, so that only their y-coordinates 
need to be tested. The children of nodes of Lf, with y- 
coordinate less than c must be considered. In particular, 
the reporting procedure proceeds recursively, as long as 
points inside the query range are found. If a point of a 
node u does not lie inside the query-range, then this point 
has y-coordinate larger than c. Therefore, all points in the 
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Fig. 1. p,: the search path, L5: the nodes that are left sons of nodes 
on P;, and do not belong to the path. 

subtree rooted at u lie outside the query-range and they 
are not reported. We can easily bound the query time 
by 0(logn + t), since O(logn) time is needed to visit 
the nodes in P^ljLi, and 0(t) time is necessary for the 
reporting procedure in their subtrees. 

E. Persistent Modified Priority Search Tree 

In this subsection, we briefly review the Modified 
Priority Search Tree of sioutas et al. [21]. Let S be a set 
of n points on the plane with coordinates (x,y), where 
x £ {1, . . . , M} and y e 5ft. Without loss of generality we 
assume that all points are distinct. We will show how to 
store the points in a data structure T, so that the t points 
in a query range of the form (— 00, b] x (—00, c], can be 
found in 0(t) time. Our structure relies on the priority 
search tree, which we augment with list-structures similar 
to those in [17]. 

We denote by T v the subtree of T with root v. 

The tree structure T has the following properties: 

> Each point of S is stored in a leaf of T and the points 
are in sorted a>order from left to right. 

> Each internal node v of T stores a point p(v) of S. 
The point p(v) is the point with the minimum y- 
coordinate amongst the points stored in the leaves 
ofT„. 

« Each node v is equipped with a secondary list S(v). 
S(v) contains the points stored in the leaves of T v 
in increasing y-coordinate. 

For convenience we will assume that the tree T is a 
complete binary tree (i.e. all its leaves have depth logn). 
Note that if n is not a power of 2, then we may add 
some dummy leaves so that T becomes complete. We 
also use an array A of size M, which stores pointers 
to the leaves of T. Specifically, A[i] contains a pointer 
to the leaf of T with maximum .x-coordinate smaller 
or equal to i. This array is used to determine in O(l) 
time the leaf of the search path P\, for b. In each leaf 
u of the tree with x-coordinate i we store the lists L(u) 
and Pl(u). The list L(u) stores the points of the nodes 
of Li. The list Pl{u) stores the points of the nodes of 
Pi which have x-coordinate smaller or equal to i. Both 
lists also contain pointers to the nodes of T that contain 



these points. Each list L(u), Pl{u), stores its nodes in 
increasing y-coordinate of their points (see Figure |2). 




Fig. 2. L(u): it stores the points of the nodes of Li. P^(u): it stores 
the points of the nodes of P, which have rr-coordinate smaller or equal 
to i. 

To answer a query of the form (— 00, b] x (— 00, c] we 
find in 0(1) time the leaf u of the search path Pf, for 
b. Then, we traverse the list Pl{u) and report its points 
until we find a point with y-coordinate greater than c. 
We traverse the list L(u) in the same manner and find 
the nodes v of whose points have y-coordinate less 
than or equal to c. For each such node v we traverse the 
secondary list S(v) and report its points until we find a 
point with y-coordinate greater than c. 

The following theorem bounds the size and the query 
time of our structure. 

Theorem 1. Given a set of n points on the plane with 
coordinates (x, y) such that x G {1, . . . , M} and y € 5ft, 
we can store them in a data structure with 0(M+n log n) 
space that supports quadrant range queries in 0(t) time, 
where t is the number of reported points. 

Proof. The query algorithm finds the t' points of 
nodes of [J Lt, that lie inside the query-range in 0(t') 
time by simple traversals of the lists Pl(u), L(u). The 
search in the respective subtrees T(v) can be performed 
by traversing the secondary lists S(v) and takes 0(t) 
additional time for reporting t points in total. Therefore, 
the query algorithm needs 0(t) time. Each list Pl{u), 
L(u) stores the respective points in the nodes of the path 
from root to u, and points in the left children of nodes 
of this path. So, the size of each list is O(logn) and the 
space of T is O(nlogn). The total space of secondary 
lists S(v) is also 0(n). Thus, the space of the whole 
structure is 0(M + nlogn) because of the size of the 
array A. 

The O(nlogn) term in the space bound is due to the 
size of the lists Pl(u) and L(u). We can reduce the total 
space of these lists to 0(n) by making them persistent. 
Ordinary structures are ephemeral in the sense that update 
operations make permanent changes to the structures. 
Therefore in ordinary structures it is impossible to access 
their old versions (before the updates). According to the 
terminology of Driscoll et al. [7] a structure is persistent, 
if it allows access to older versions of the structure. There 
are two types of persistent data structures: partially and 
fully persistent. A partially persistent data structure allows 



updates of its latest version only, while a fully persistent 
one allows updates of any of its versions. In [7], a general 
methodology is proposed for making data structures of 
bounded in-degree persistent. With their method such 
a structure can be made partially persistent with 0(1) 
amortized space cost per change in the structure. In our 
case a list can be made partially persistent with a 0(1) 
amortized increase in space per insertion/deletion. 

We show how to implement the lists Pl(u) using a 
partially persistent list. Let u be a leaf in T and let w be 
its predecessor (the leaf on the left of u). We denote by 
x u the ^-coordinate of u and by x w , the ^-coordinate 
of w. The two root-to-leaf paths P Xu , P Xm , share the 
nodes from the root of T to the nearest common ancestor 
of u, w. As a result, we can create Pl{u) by updating 
Pl(w) in the following way. First we delete from list 
Pl{w) the points that don't lie on P Xu . Then we insert the 
points of P Xu which have x-coordinate smaller or equal 
to x u . In this way we can construct all lists as versions 
of a persistent list: we begin from the leftmost leaf and 
construct the list Pl(u) of each leaf u by updating the 
one of its predecessor (see Figure |3). 




Fig. 3. Lists Pl{u) and L(u) are implemented as partially persistent 
lists, by performing a sweep from left to right. 

The total number of insertions and deletions is 0{n) 
because each point is inserted and deleted only once. 
Therefore the space of all the lists is 0(n). In the same 
way, lists L(u) are constructed for all leaves in 0(n) 
space. Therefore: 

Theorem 2. Given a set of n points on the plane with 
coordinates in the range [1, M] x R we can store them in a 
data structure with 0(n + M) space that allows quadrant 
range queries to be answered in 0(t) time, where t is the 
number of answers. 

The preprocessing time is 0(M + nlogn) but with a 
more careful implementation we can reduce this complex- 
ity to 0(M + n), by using the pruning technique as in 
[12]. 

F. Geometric Transformation Of Duality 

In homogeneous coordinates, there exists a duality 
between the point (a, b, c) and the line (ax + by + cz = 0) 
for all (a, b, c) ^ (0,0,0). Furthermore, the following 
pairs are dual: 

• points on a line < — ► lines throught a point 



< line segment < — ► "double wedge" of lines 

< set of lines intersecting a line segment < — ► set of 
points in a double wedge (see Figure |4]l 

A more formal definition of this powerfull tool can be 
found in [5]. 




Fig. 4. An example of Geometric Duality Transformation: Dual pairs 

III. Algorithms for Orthogonal Objects 

A. Data Structures for Orthogonal Triangle Range 
Queries 

An orthogonal triangle range query is the problem of 
determining all the points from a set S of n points on 
the plane lying inside an orthogonal triangle. Recall, a 
triangle is orthogonal iff two of its edges are axis-parallel. 
Let T be an orthogonal triangle defined by the point 
(x q ,y q ) and the line L q that is not axis-parallel (see 
Figure |5j. 




Fig. 5. An example of Orthogonal Triangle Range Query 

A retrieval query for this problem can be supported 
efficiently by the following data structures: 

1) First Solution: A 3-layered Tree: To set up the 

data structure, first sort the n points according 
to their .x-coordinates and then store the ordered 
sequence in a leaf-oriented balanced binary search 
tree of depth O(logn). This structure answers the 
query: "determine the points having x-coordinates 
in the range [ij,^]" by traversing the two paths 
to the leaves corresponding to x\, X2. The points 
stored as leaves at the subtrees of the nodes which 
lie between the two paths are exactly these points 
in the range [xi,X2]- For each subtree, the points 
stored at its leaves are organized further to a second 
level data structure according to their y-coordinates 
in the same way. For each subtree of the second 
level data structure, the points stored at its leaves 
are organized further to a third level modified data 
structure of Chazelle (presented in subsection 2.3) 



for half-plane range query. So, each Orthogonal 
Triangle Range Query is performed through the 
following steps: 

a) In the tree storing the pointset S according 
to ir-coordinates, traverse the path to x q . All 
the points having ^-coordinate in the range 
[x q , oo ) are stored at the subtrees on the nodes 
that are right sons of a node of the search path 
and do not belong to the path. There are at 
most [logn] such disjoint subtrees. 

b) For every such subtree traverse the path to y q . 
By a similar argument as in the previous step, 
at most [logn] disjoint subtrees are located, 
storing points that have y-coordinate in the 
range [y q ,oo). 

c) For each subtree in Step 2, apply the half- 
plane range query algorithm presented in sub- 
section 2.3 in order to retrieve the points that 
lie on the side of line L q towards the triangle. 

The correctness of the above algorithm follows from 
the data structure used. Since we have to visit 
O(logn) subtrees in each of the two first layers 
and the third layer requires 0( l ^°^ 1 +A) time, the 

overall query time becomes 0( l ° g 9 lo ^ n + A) while 
the space complexity is 0(nlog 2 n). 
2) Second Solution: A two layered tree: The first 
layer is a linear space Persistent Modified Priority 
Search Tree. The second layer organizes each of 
the 0(log 2 n) subsets of Pl(u) and L(u) with the 
modified structure of [5] presented in subsection 2.3 
which returns a set of nodes v whose points lie 
inside the triangle. Then, it organizes each of the 
0(n 2 ) subsets of S(v) with the same structure (see 
figure 6). 

More specifically, we find in O(l) time the leaf 
u of the search path Pb for b. Each leaf u stores 
two y-ordered Persistent lists Pl{u) and L{u), both 
of size 0(logn). The total number of subsets, 
which can be produced by taking all the possible y- 
ordered permutations amongst the y-coordinates of 
Pl(u) and L(u), is O(logn), since each list Pl(u) 
or L(u) consumes O(logn) space in worst-case. 
Each subset has minimum size 1 and maximum 
0(log n), thus their total space is 1 + 2 + + logn = 
0(log 2 n). We can access the appropriate subset of 
points (according to y q - coordinate of the query) 
in O (log logn) time by applying a simple binary 
searching. Then we apply in the subset above a half 
plane range query (according to query line L q ) in 
order to find the nodes v whose points lie inside 
the triangle. For each such node v we organize its 
secondary list S(v) again as the modified struc- 
ture of Chazelle presented in subsection 2.3. The 
total number of subsets, which can be produced 
by taking all the possible y-ordered permutations 
amongst the y-coordinates of S(v), is 0(n), since 
each secondary list S(v) consumes 0(n) space in 



worst-case. Each subset has minimum size 1 and 
maximum 0(n), thus their total space is 1 + 2 + + 
n = 0(n 2 ). We can access the appropriate subset 
of points (according to y q - coordinate of the query) 
in 0( , lo P n ) time by using the fusion tree method 

v loglogn ' J id 

[11]. Then we apply in the subset above a half 
plane range query (according to query line L q ). 
As a consequence the overall query time becomes 
0(logn/ loglogn + A) and the space consumption 
is 0(n 2 ). 




Fig. 6. The second solution. The lists P, L, S and its respective Half 
Plane (HP) structures 



B. The Results 

Every canonical fc-vertex polygon can be decomposed 
into 0(k) orthogonal triangles because of the strict sym- 
metry in the topology of the vertices. In Figure [7] such a 
decomposition is depicted. 
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Fig. 7. Orthogonal triangle decomposition of a canonical octagon 

• Solution 1: Querying the 3-layered data struc- 
ture 0(k) times, the 0(n log 2 n) space and 
0(klog 3 n/loglogn + A) time follows. 

• Solution 2: Querying the two layered data 
structure 0(k) times, the 0(n 2 ) space and 
Oiklognj loglogn + A) time follows. 




IV. Some Thoughts on the General Polygon 
Retrieval Problem 

A. The Algorithm of Paterson and Yao 

This section sketches the basic ideas of the algorithm 
for polygon retrieval on points that Paterson and Yao [13] 
have presented. 

This problem is to preprocess a set S of n points on 
the plane, so that for any query polygon P the subset 
of them lying inside it can be reported efficiently. We 
assume that the given polygon is convex with fc vertices. If 
not, (P is non-convex) then a suitable decomposition into 
convex parts can be carried out and the algorithm can be 
applied to every such part. The key-idea of Paterson/Yao's 
algorithm is a further decomposition of every convex part 
into O(fc) simpler parts called quads. They solve the 
problem for each such quad separately, using the well- 
known geometric transformation of duality. 

To set up the data structure, Paterson/Yao first sort the 
n points according to their polar angles at the origin and 
then store the ordered sequence in a leaf-oriented bal- 
anced binary search tree of depth 0(log n). This structure 
answers the query: "determine the points having polar 
angle in the range [a\ , (12] by traversing the two paths to 
the leaves corresponding to ax, 0,2- The points stored as 
leaves at the subtrees of the nodes which lie between the 
two paths are exactly these points in the range [ax, 0,2]. For 
each subtree, the points stored at its leaves are organized 
further to a second level data structure as follows: by the 
duality correspondence, these points are mapped to a set 
of straight lines, and a point location structure is built 
for this planar subdivision [8]. Using this data structure 
a query of the form: "Report the points lying in a double 
wedge" is reduced to a query to the dual structure of the 
form: "report the lines that intersect a query line segment." 

At this point, we are ready to give a step by step 
description of the polygon retrieval algorithm of Paterson 
and Yao: 

1) Let P be the fc-vertex convex query polygon. Sepa- 
rate P into simpler regions called quads by cutting 
the plane into sectors. The sectors are obtained by 
drawing semi-infinite lines from the origin to each 
vertex of P. Each quad has a particularly simple 
form (it is the intersection of a sector and a double 
wedge), (see Figure [8) 




Fig. 8. An example of Sectors and Quads 

Every fc-gon is the disjoint union of O(k) such 
quads. 

2) Using the two-level data structure do the following 
for each quad with bounding rays l\, £2' 



• Determine the (at most) 2 \logn\ disjoint sub- 
trees that descend from the two search paths 
to the angles of £\, £2 and lie between them. 
All the points between l\ and I2 are stored as 
leaves in these subtrees. 

• for each such subtree in turn report the points 
inside the double wedge determining the quad 
using the dual point location structure. 

The space used is 0(n 2 ) because of the point location 
structure and the query time is 0(k\ogn + A) (there 
exist O(fc) quads and the time spent for every "quad- 
retrieval" is 0(log 2 n + Ai) which can be improved 
using the fractional cascading technique [4] on the point 
location structures). It is worthwile to note that the above 
algorithm works even when P is unbounded. If P is 
non-convex, then the preliminary decomposition step into 
convex parts does not change the asymptotic space and 
time bounds. 

Unfortunately, the 0(n 2 ) space bound is the best pos- 
sible if we directly use the dual point location structures. 
This means that in order to reduce the space requirements 
of the later algorithm, we have to find an algorithm 
supporting "quad retrieval" without using the duality 
on arbitrary sets of 0(n) points which lies to 0(n 2 ) 
space subdivisions. The next section demonstrates such 
an algorithm working efficiently when P is a canonical 
fc-vertex polygon. 

B. An Extended Approach 

In this section we present an extended approach for the 
general polygon retrieval problem. The proposed solution 
is based on the algorithm of Paterson/Yao for answering 
"quad retrieval" queries, which was described in the 
previous section. The key-idea is the reduction of a quad 
retrieval query to a constant 0(1) number of orthogonal 
triangles retrieval queries. This is achieved by drawing 
parallel lines (vertical or horizontal) from the vertices 
of the quad so that the lines stay inside the quad and 
intersect the rays from the origin. To do this, a suitable 
decomposition of each quad is needed. In the general case 
of an arbitrary (non canonical) fc-vertex polygon, such 
a decomposition is not always possible. For example, 
the quad in Figure |9ja) can be decomposed into 0(1) 
orthogonal triangles while the quad in Figure |9jb) cannot 
be decomposed. 

Not part of 




(a) (b) 



Fig. 9. An example of Quad Decomposition 

This of course depends on various geometric character- 
istics of the quad. If the polygon is fc-vertex canonical, we 



can always directly carried out a decomposition into O(k) 
orthogonal triangles because of the strict symmetry in 
the topology of the vertices. Arbitrary fc-vertex polygons 
satisfy such a decomposition only in special cases. 

V. Conclusions 

In this work we presented efficient algorithms for the 
problem of canonical polygon retrieval queries on the 
plane. One deficiency or our algorithms is that they 
are static. It seems very interesting the dynamization 
of presented algorithms. Furthermore, it still remains 
an open problem whether an 0(nlog°^ra) space and 
0(log°^ n + A) time algorithm exists for the general 
problem of arbitrary polygon retrieval queries. Finally, 
it would be of great importance the incorporation of 
the proposed algorithms in meandric polygons [19], a 
specific kind of planar polygons with many applications 
in bioinformatics. The latter constitutes a real case study 
in bioinformatics, which we plan to implement in the near 
future. 
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